CRDKNSWC-SSM-64 — 92/22  Measurement  of  the  Out-of-Plane  Shear  Response  of  Thick  Section 

Composite  Materials  Using  the  V-notched  Beam  Specimen 


AD-A264 


Carderock  Division 

Naval  Surface  Warfare  Center 


Bethesda.  Md.  20084-^0 


CRDKNSWC-SSM-64— 92/22  April  1993 

Survivability,  Structures,  and  Materials  Directorate 
Technical  Report 


Measurement  of  the  Out-of-Piane  Shear  Response 
of  Thick  Section  Composite  Materials  Using  the 
V-notched  Beam  Specimen 


by 

K.  Gipple 
D,  Hoyns 


SDTIC 

ELECTE  1^ 
MAY  2  4  1993  I  I 

A  ^ 


93-11445 


Carderock  Division 

Naval  Surface  Warfare  Center 

Bethesda,  Md.  20084-5000 


CRDKNSWC-SSM-64— 92/22  April  1993 

Survivability,  Structures,  and  Materials  Directorate 

Technical  Report 


Measurement  of  the  Out-of-Plane  Shear  Response 
of  Thick  Section  Composite  Materials  Using  the 

V-notched  Beam  Specimen 

by 

K.  Gipple 

■  ^  D.  Hoyns 


Approved  for  piArlic  release;  distrfl>utlon  Is  unlimited. 


TABLE  OF  CONTENTS 


List  of  Tables .  iii 

List  of  Figures .  iv 

Abstract .  1 

Administrative  Information .  1 

IntroducdcHi .  1 

Background:  Moire  Interferometry .  2 

Experimental .  3 

Materials  and  Specimen  Preparation .  3 

Instrumentation  and  Mechanical  Testing .  4 

Conventional  Method .  4 

Moird  Interferometry  and  New  Full  Section  Gage .  4 

Microscopy .  5 

Three  Dimensional  Hnite  Element  Analysis .  5 

Finite  Element  Analysis  Program  Description .  5 

Nonlinear  Material  Model .  5 

Finite  Element  Model  of  the  V-notched  Beam  Specimen .  7 

Finite  Element  Computations .  8 

Results .  9 

Experimental  Results .  9 

Conventional  Gages .  9 

Moir6 .  9 

Moire/Full  Section  and  Back-to-Back  Full  Section  Gages .  1  0 

SEM .  1  1 

Experimental/Analytical  Results .  1  1 

Analysis  Versus  Conventional  Gage  Results .  1  1 

Analysis  Versus  Moire/Full  Section  Gage  Data .  1  2 

Conclusions .  1  2 

Acknowledgments .  1  3 

Refo^nces .  1  4 


LIST  OF  TABLES 


Page 

1.  Out-of-plane  shear  modulus  and  strength  values  for 

unidirectional  and  crossply  AS4/3501-6  and  S2  glass/3501 -6 


laminates . .  17 

2.  AS4  graphite  and  S2  glass  Hber  properties  used  in  analysis .  18 

3.  3501-6  epoxy  matrix  properties  used  in  simplified  unit  cell 

model .  19 


LIST  OF  FIGURES 


Page 


1.  Schematic  of  moir6  setup  after  Guo  and  Ifju  [15] .  20 

2.  V-notched  beam  specimen  (a)  specimen  geometry  and 

(b)  Hnite  element  mesh  and  boundary  conditions .  21 


3.  Strain  gages  used  in  study  (a)  conventional  0/90  stacked  gage, 

(b)  full  section  gages  in  st^^tked  and  side-by*side  configurations....  22 


4.  V-notched  beam  test  specimen  mounted  in  test  fixture  modeled 

after  fixture  developed  by  Adams  and  Walrath  [11] .  23 

5.  Simplified  unit  cell  model  based  on  spring  analogy  showing  the 

load  path  variations  between  fiber  and  matrix  dominant 
components . .  24 

6.  Predicted  through-thickness  shear  behavior  for  unidirectional 
(G23)  AS4/3501-6  and  S2  glass/350i-6  using  WY02D  and 

the  simplified  unit  cell  model .  25 

7.  Predicted  through-thickness  shear  behavior  for  AS4/3501-6  and 
S2  glass/3501-6  crossply  laminates  using  the  simplified  unit 

cell  model .  25 

8.  U  (a)  and  V  (b)  moir€  field  fringe  patterns  for  unidirectional 

(G23)  AS4/3501-6  specimens  loaded  at  23  MPa  and  7.7  MPa 
respectively .  26 

9.  U  (a)  and  V  (b)  moire  field  fringe  patterns  for  AS4/3501-6 

crossply  specimen  at  a  load  level  of  14.3  MPa .  27 

10.  U  (a)  and  V  (b)  moire  field  fringe  patterns  for  unidirectional 

(G23)  S2  glass/3501-6  specimens  loaded  at  21  MPa .  28 

11.  U  (a)  and  V  (b)  moir6  field  fringe  patterns  for  S2  glass/3501-6 

crossply  specimens  loaded  at  21.8  h^a  and  11.3  MPa 
respectively .  29 


12.  Normalized  shear  strain  distribution  along  the  test  section  of  a 
unidirectional  (G23)  AS4/3501-6  V-notched  beam  specimen 


IV 


measured  by  high  sensitivity  moirg  at  two  load  levels .  30 

13.  Normalized  shear  strain  distribution  along  the  test  section  of  a 

unidirectional  (G23)  S2  glass/3501-6  V-notched  beam  specimen 
measured  by  high  sensitivity  moire  at  several  load  levels .  30 

14.  Normalized  shear  strain  distribution  along  the  test  section  of  a 
crossply  AS4/3501-6  V-notched  beam  specimen  measured 

by  high  sensitivity  moir6  at  two  load  levels .  31 

15.  Normalized  shear  strain  distribution  along  the  test  section  of  a 
crossply  S2  glass/3501-6  V-notched  beam  specimen  measured 

by  high  sensitivity  moire  at  two  load  levels .  31 

16.  Shear  stress-strain  response  for  AS4/3501-6  and 

S2  glass/3501-6  unidirectional  (G23)  specimens  instrumemed 

with  conventional  gages  vs.  moire  and  full  section  gage  data .  32 

17.  Shear  stress-strain  response  for  AS4/3501-6  and 

S2  glass/3501-6  crossply  specimens  instrumented  with 
conventional  gages  vs.  moir6  and  full  section  gage  data .  32 

18.  Photomicrographs  of  AS4/3501-6  and  S2  glass/3501-6 

crossply  specimen  failures . 33 

19.  Photographs  of  AS4/3501-6  and  S2  glass/3501-6 

unidirectional  failed  specimens .  34 

20.  Shear  stress-strain  response  for  unidirectional  (G23) 

AS4/3501-6  and  S2  glass/3501-6  specimens  instrumented 
with  conventional  gages  (062)  vs.  predicted  values  based  on 
finite  element  analyses  with  linear  and  nonlinear  elastic 

material  properties . 35 

21.  Shear  stress-strain  response  of  crossply  specimens  as 

predicted  by  finite  element  analyses  with  linear  and 
nonlinear  elastic  material  properties,  and  measured  by 
conventional  (062)  gages .  35 

22.  Average  shear  stress-strain  response  for  AS4/3501-6  and 
S2  glass/3501-6  unidirectional  (G23)  specimens  measured 
with  moird  interferometry  and  full  section  gages  vs.  predicted 
values  based  on  finite  element  analyses  with  nonlinear  elastic 


v 


material  properties. 


36 


23,  Average  shear  stress-strain  response  for  AS4/3501-6  and 
S2  glass/3501-6  crossply  specimens  measured  with  moir6 
interferometry  and  full  section  gages  vs.  predicted  values 
based  on  finite  element  analyses  with  nonlinear  elastic 
material  properties . 


36 


VI 


ABSTRACT 


The  out-of-plane  shear  response  of  thick,  unidirectional 
and  crossply,  AS4I3S01-6  and  S2  glass/ 3 50 J -6  laminates  was 
investigated  theoretically  and  experimentally.  The  specimens 
considered  were  of  the  V-notched  beam  (losipescu) 
configuration.  Strains  were  monitored  in  the  specimen  test 
sections  using  conventional  strain  gages,  moir6 
interferometry  and  new  full  section  strain  gages.  Crossply 
laminates  exhibited  a  fairly  uniform  strain  distribution  away 
from  the  notches  resulting  in  agreement  of  initial  moduli 
obtained  through  all  three  strain  measurement  methods. 
Strains  measured  with  conventional  gages  centered  between 
he  notches  of  the  unidirectional  specimens  exceeded  the 
average  strain  from  notch  to  notch,  as  measured  by  moiri 
interferometry  and  the  full  section  gages.  Therefore,  a  lower 
initial  modulus  was  calculated  using  conventional  gage  data 
vs.  moiri  or  full  section  data.  Measured  strains  correlated  very 
well  with  predicted  strains  from  specimen  finite  element 
analyses  using  nonlinear  elastic  material  properties. 
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INTRODUCTION 

Increasing  the  efficiency  of  composite  designs  requires  a  better 
understanding  of  the  three-dimensional  (3-D)  behavior  of  composite 
materials.  Of  particular  interest  is  the  oul-of-plane  or  throu  nickness 
direction.  The  combination  of  through-thickness  shear  and  sile  stresses 
with  the  inherently  low  through-thickness  strength  of  composite  mat'^rials 
plays  a  vital  role  in  composite  structural  behavior  [1-6].  At  some 
geometric  discontinuities  such  as  free  edges  or  ply  dropoff  regions, 
through-thickness  stresses  can  be  minimized  by  careful  choice  of  lamina 
stacking  sequence  [7]  or  proper  design  procedures  [3].  However,  in  many 
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instances,  the  significance  of  through -thickness  stresses  must  be  assessed 
on  a  case  by  case  basis  due  to  fixed  laminate  or  geometry  requirements  [3]. 
Three  dimensional  analyses  of  structures  or  structural  details  dominated 
by  through-thickness  stresses  are  also  complicated  by  the  scarcity  of 
through-thickness  material  data  due  to  the  lack  of  suitable  specimens  tc 
measure  out-of-plane  behavior  [5,8]. 

The  current  study  focuses  on  the  measurement  of  the  through¬ 
thickness  shear  response  using  the  V-notched  beam  (losipescu)  specimen 
in  the  modified  Wyoming  test  fixture  [9].  The  V-notched  beam  specimen  is 
typically  used  for  in-plane  shear  testing  of  composites.  It  has  been 
rigorously  investigated  with  regards  to  specimen  geometry  and  loading 
effects  [10],  plasticity  and  material  orientation  effects  [11],  strain  field 
uniformity  and  purity  [10],  and  repeatability  of  results  [12].  Adams  et  al. 
[13]  did  use  the  V-notched  beam  specimen  for  through-thickness  shear 
measurements  comparing  the  behavior  of  unidirectional  and  cross-ply 
laminates  fabricated  by  bonding  several  thin  laminates  together  vs.  those 
cut  from  a  single  thick  laminate.  To  the  authors'  knowledge,  no  study  has 
investigated  the  use  of  the  losipescu  specimen  for  through-thickness  shear 
testing  more  throroughly. 

In  this  program,  moire  interferometry  and  3-D  finite  element  analyses 
(3-D  elastic  and  3-D  nonlinear  elastic  material  properties)  are  used  to 
evaluate  the  uniformity  and  purity  of  the  strain  field  in  the  gage  section  of 
the  through-thickness  shear  V-notched  beam  specimen.  The  impact  of 
nonuniform  fields  on  strain  gage  measurements  is  also  discussed.  The 
results  from  an  experimental  gage  developed  at  the  Virginia  Polytechnic 
Institute  and  State  University  (VPI&SU)  [14]  are  included  for  comparison 
to  conventional  gages  and  to  assess  the  differences  in  average  strain 
readings  between  the  front  and  the  back  of  the  specimen. 

BACKGROUND:  MOIRE  INTERFEROMETRY 

High  sensitivity  moir6  interferometry  is  an  optical  technique  for  whole 
field  deformation  measurement.  In  this  study,  a  collimated  beam  of  light 
was  directed  to  a  three  mirror  interferometer  (Figure  1)  producing  two 
orthogonal  virtual  reference  gratings  of  frequency  f  [15].  These  gratings 
interact  with  a  cross-lined  diffraction  grating  of  frequency  f/2,  replicated 
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onto  the  specimen  surface,  to  create  contour  maps  of  U  (horizontal)  and  V 
(vertical)  displacement  fields.  These  contour  maps  are  composed  of 
alternating  black  and  white  fringe  patterns.  The  surface  normal  strains,  Cx 
and  Ey,  and  shear  strain,  Yxy.  are  related  to  the  displacement  fields  by  the 
equations: 

ax  flax]  (1) 


_3V_l(3Ny\ 

y  By  f\  dy  / 

(2) 

Yxy 


_  au  ^  av  _  ij 

By  dx  f 


dy  dx 


(3) 


where  Nx  and  Ny  correspond  to  the  U  and  V  field  fringe  order, 
respectively. 


EXPERIMENTAL 
Materials  and  Specimen  Preparation 

The  material  systems  evaluated  in  this  investigation  were  Hercules 
AS4/3501-6  graphite/epoxy  and  S2/3501-6  glass/epoxy  prepreg.  Thick- 
section  panels  of  unidirectional  [0]i40  and  crossply  10/90]72S  construction 
were  fabricated  and  autoclave  cured  at  CARDEROCKDIV,  NSWC.  As 
specified  in  Figure  2,  V-notched  beam  specimens  were  machined  from  the 
unidirectional  panels  in  the  90°  direction,  and  will  be  referred  to  as  the 
G23  specimens.  Specimens  from  the  crossply  panel  were  machined  in  the 
QO  and  90®  directions  to  obtain  Gxz  and  Gyz-  Unidirectional  AS4  specimens 
were  0.51  cm.  (0.20  in.)  thick,  unidirectional  S2  glass  specimens  were  0.38 
cm.  (0.15  in.)  thick,  and  crossply  specimens  of  both  materials  were  0.53 
cm.  (0.21  in.)  thick.  Fiber  volumes  were  determined  by  standard  matrix 
digestion  and  incineration  methods,  ASTM  D3171  [16]  and  ASTM  D2584 
[17]  respectively.  Specimen  surfaces  were  ground  and  polished  with  1200 
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micron  grit  paper  using  standard  metallographic  techniques. 

Instrumentation  and  Mechanical  Testing 
Co"  entional  Method 

Stacked  (0O/90O)  strain  rosettes  (Measurements  Group  CEA-06- 
062WT-350)  were  centered  between  the  notch  roots  at  +45/-45  degrees 
to  the  loading  axis  (Figure  3).  Gages  were  bonded  to  both  the  front  and 
back  faces  of  the  specimens.  Compliant  shims  were  also  placed  along  the 
edges  of  the  specimens  to  avoid  possible  twisting  from  localized 
asperities,  as  reported  by  Morton,  et  ai,  [10].  The  CARDEROCKDIV,  NSWC 
V-notched  beam  test  fixture,  based  on  specifications  for  the  modified 
Wyoming  version  developed  by  Walrath  and  Adams  [11  ]»  is  pictured  in 
Figure  4.  Specimens  were  centered  in  the  fixture  and  loaded  to  failure 
at  1.25  mm/min.(,05  in./min.)  by  a  hydraulic,  267  kN  (60,000  lb.) 
capacity,  Satec  testing  machine.  Strain  output  was  conditioned  through 
2310  amplifiers  and  recorded  by  a  microcomputer  based  data 
acquisition  system. 

Moir6  Interferometry  and  New  .  Full  Section  Strain  Gage 

Specimen  gratings  approximately  0.0254  mm  (0.001  in.)  thick  were 
replicated  onto  the  front  surface  of  specimens  of  each  material  using  the 
method  described  in  Reference  18.  The  previously  described  moire 
interferometry  system  was  set  up  on  the  cross-head  of  the  test  machine 
and  produced  a  virtual  grating  frequency  of  2400  1/mm  (60,960  1/in). 
The  corresponding  sensitivity  of  the  displacement  measurements  is  1/f, 
or  0.417p,m(16.4pin.)  per  fringe  order.  The  application  of  gratings  and 
the  interferometry  testing  were  performed  at  VPI&SU. 

New  full  section  strain  gages  with  side  by  side  ■♦-45®/-450 
configuration  were  developed  at  VPI&SU  and  recently  made 
commercially  available.  One  of  these  gages  was  bonded  to  the  back 
surface,  opposite  the  moir6  grating,  of  each  of  the  four  specimens. 
Schematics  of  the  gages  used  in  this  study  are  shown  in  Figure  3.  The 
loading  fixture  was  of  the  modified  Wyoming  type,  and  the  testing 
machine  was  a  screw-driven  267  kN  (60,000  lb.)  capacity  Tinius  Olsen. 
U  and  V  displacement  fields  in  the  test  section  of  the  specimens  were 


4 


recorded  on  pan  film  at  several  load  levels,  from  the  clamped  no  load 
condition  to  failure. 

Two  additional  unidirectional  (G23)  and  crossply  specimens  of  each 
material  were  also  instrumented  with  back-to-back  full  section  gages. 
These  specimens  were  tested  with  the  CARDEROCKDIV,  NSWC  fixture  and 
testing  machine  referenced  previously.  Table  1  shows  the  number  of 
specimens  tested  for  each  orientation,  material  and  strain  measurement 
method. 

Microscopy 

Failed  crossply  specimens  were  sputter  coated  with  Au/Pd,  then 
inspected  via  a  JEOL  model  35CF  scanning  electron  microscope. 
Photomicrographs  were  taken  of  failure  surfaces. 

THREE-DIMENSIONAL  FINITE  ELEMENT  ANALYSIS 
Finite  Element  Analysis  Program  Description 

A  3-D  finite  element  analysis  of  the  V-notched  beam  through-thickness 
shear  coupon  was  conducted  using  ABAQUS  [19].  ABAQUS  is  a  general 
purpose  finite  element  program  developed  by  Hibbitt,  Karlsson,  and 
Sorensen,  Inc.  It  provides  a  wide  range  of  options  for  numerical  modeling 
of  structural  response.  Stress  analyses  may  be  static  or  dynamic,  linear  or 
nonlinear.  Sources  of  nonlinearity  include  material,  geometric  and 
boundary  condition  nonlinearities.  A  number  of  nonlinear  material  models 
are  available  in  ABAQUS  or  the  material  may  be  described  by  the  user  in  a 
UMAT  subroutine.  The  latter  option  was  used  in  this  analysis  which 
considered  the  effects  of  both  linear  and  nonlinear  elastic  material 
behavior  on  the  response  of  glass  and  graphiie/epoxy,  unidirectional  and 
crossply  V-notched  beam  through -thickness  shear  specimens. 

Nonlinear  Material  Model 

Development  of  the  nonlinear  material  model  used  in  the  analysis 
followed  the  framework  outlined  in  Reference  20.  This  framework  consists 
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of  three  steps:  development  of  the  detailed  3-D  nonlinear  material 
response  using  micromechanical  finite  element  analyses,  adjustment  of  a 
simplified  unit  cell  material  model  to  match  the  detailed  analysis  results, 
and  transition  of  the  simplified  model  to  a  general  purpose  finite  element 
program  for  structural  analyses. 

Detailed  finite  element  analyses  of  a  fiber-matrix  unit  cell  were 
performed  using  WY02D  [21]  to  characterize  the  nonlinear  composite 
stress-strain  behavior.  WY02D  is  a  2-D,  generalized  plane  strain  finite 
element  analysis  program  developed  by  the  Composite  Materials  Research 
Group  in  the  Department  of  Mechanical  Engineering  at  the  University  of 
Wyoming.  It  is  specifically  tailored  for  micromechanical  analyses  of 
continuous  fiber-reinforced  composite  materials. 

Transversely  isotropic  fiber  properties  and  elastoplastic  matrix 
properties  were  used  in  the  micromechanical  analyses.  The  plastic 
response  of  the  matrix  is  initiated  with  an  octahedral  shear  stress  yield 
criterion,  and  governed  by  the  Prandtl-Reuss  flow  rule  which  assumes  that 
plastic  strain  is  proportional  to  the  deviatoric  stress  [22,23].  Matrix 
octahedral  shear  stress-strain  behavior  is  developed  from  Richard- 
Blacklock  cu've  fits  of  matrix  solid  rod  torsion  data  at  various 
temperatures  and  moisture  levels, 

A  simplified  unit  cell  model,  developed  by  Pecknold  [24],  was  chosen  to 
approximate  the  micromechanical  finite  element  analysis  results  and  to 
interface  to  ABAQUS.  The  unit  cell  configuration  is  based  on  fiber  and 
matrix  "spring”  elements  in  series  and/or  in  parallel  depending  on  the  load 
path  through  these  elements  (see  Figure  5).  As  in  WY02D,  the  fiber  is 
considered  transversely  isotropic.  In  contrast  to  the  elastoplastic  matrix 
material  model  in  WY02D,  the  matrix  shear  behavior  is  nonlinear  elastic 
and  expressed  in  a  Ramberg-Osgood  form  [25].  The  fiber  and  matrix 
properties  used  in  this  model  are  listed  in  Tables  2  and  3.  Stress  averaging 
techniques  are  used  to  develop  the  response  of  the  unidirectional  lamina 
from  the  fiber-matrix  unit  cell.  Unidirectional  lamina  are  assembled  into 
sublaminates  (smallest  repeating  unit  within  the  laminate)  based  on 
classical  lamination  theory  and  the  condition  that  out-of-plane  stresses  at 
the  layer  interface  must  be  continuous  across  the  interface  to  satisfy 
equilibrium.  The  inhomogenous  sublaminate  is  replaced  by  a 
homogeneous  anisotropic  material  using  equivalent  continuum  modeling 
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techniques. 

Initial  comparisons  of  the  predicted  in-plane  shear  response  using  the 
simplified  unit  cell  model  and  WY02D  with  matching  matrix  material 
models  based  on  neat  resin  data  showed  that  the  simplified  unit  cell  model 
predicted  a  stiffer  response.  The  matrix  material  parameters  in  the 
simplified  unit  cell  model  were  adjusted  so  that  its  predicted  in-plane 
shear  response  matched  WY02D‘s  for  both  AS4/3501-6  and  S2 
glass/3501-6  [20].  Figure  6  compares  the  predicted  unidirectional 
composite  through-thickness  shear  stress-strain  behavior  using  WY02D 
and  the  simplified  unit  cell  model  with  the  adjusted  matrix  material  for 
AS4/3501-6  and  S2glass/3501-6.  There  is  good  agreement  between  both 
predicted  responses  for  AS4/3501-6,  but  the  simplified  unit  cell  model 
predicts  a  softer  through-thickness  shear  response  for  S2  glass/3501 -6. 
This  effect  has  not  been  investigated,  but  may  be  related  to  the  use  of 
average  rather  than  local  stresses  in  the  simplified  unit  cell  model  [24]. 
Figure  7  shows  the  predicted  behavior  for  AS4/3501-6  and  S2  glass/3501 - 
6  crossply  laminates  using  the  simplified  unit  cell  material  model  only.  A 
composite  laminate  analysis  is  not  possible  with  the  micromechanical  unit 
cell  in  WY02D  since  it  can  be  used  to  represent  unidirectional  composite 
geometries  only.  Since  these  curves  represent  the  material  response  to  a 
pure  and  uniform  through-thickness  shear  loading,  they  are  not 
necessarily  expected  to  match  the  response  generated  with  the  V-notchel 
beam  specimen.  The  predicted  material  response  is  used  in  a  finite 
element  analysis  o^  the  specimen  with  its  particular  boundary  and  loading 
conditions,  to  generate  the  predicted  response  for  the  specimen. 

Finite  Element  Model  of  the  V-notched  Beam  Specimen 

A  3-D  finite  element  model  was  constructed  following  the  dimensions 
of  the  V-notched  beam  specimen  in  Figure  2.  Eight  noded  solid  continuum 
brick  elements  (C3D8  in  ABAQUS  notation)  were  used  throughout  the 
model.  The  total  number  of  elements  was  250  with  one  element  through 
the  thickness.  Maximum  element  aspect  ratios  were  3.6  or  less  for 
lengthiwidth  and  6.25  or  less  for  length: thickness. 

Figure  2  also  illustrates  the  boundary  conditions  applied  to  the  model. 
On  one  side  of  the  specimen,  the  nodes  along  the  specimen  edge  in  contact 
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with  the  fixture  were  fixed  in  the  vertical  direction.  A  vertical 
displacement  was  imposed  on  the  other  side.  These  are  highly  idealized 
boundary  conditions.  The  effects  of  loading  variations  on  solution  results 
has  been  addressed  for  in<plane  shear  specimens  by  Morton  [10],  but  will 
not  be  considered  here. 

Finite  Element  Computations 

The  simplified  material  model  is  interfaced  to  ABAQUS  through  a 
UMAT  subroutine  written  by  Hajali  [26].  Strains  are  updated  directly 
within  the  algorithm  by  adding  the  strain  increment  to  the  total  strain  at 
the  previous  increment.  Both  the  total  strain  at  the  previous  increment 
and  the  strain  increment  for  the  equivalent  sublaminate  material  are 
supplied  to  the  subroutine  by  ABAQUS.  Micro  level  stress  and  strain 
increments  are  approximated  using  sublaminate  stress  and  strain 
increments  and  tangent  material  properties  in  a  reverse  application  of 
Pecknold's  unit  cell  to  sublaminate  model  described  in  the  preceding 
section.  The  initial  sublaminate  stress  increment  will  satisfy  equilibrium 
and  compatability,  but  not  the  nonlinear  constitutive  equations  at  the 
micro  level.  The  stress  increment  is  adjusted  in  an  iterative  process 
(modified  Newton-Raphson  method)  until  the  constitutive  conditions  are 
satisfied  to  a  specified  tolerance.  Future  refinements  of  the  stress  update 
procedure  are  expected  to  improve  computational  speed  and  convergence 
time 

In  ABAQUS,  a  modified  Riks  algorithm  was  chosen  to  define  the  load  on 
the  specimen  throughout  the  load  step.  This  algorithm  assumes 
proportional  loading  and  that  the  load  magnitude  in  each  increment  of  the 
step  is  part  of  the  solution  rather  than  being  specified  by  the  user. 
Convergence  within  each  increment  was  based  on  an  acceptable  tolerance 
measure  for  the  solution  of  the  equilibrium  equations.  ABAQUS 
recommends  a  tolerance  measure  two  to  four  magnitudes  less  than  typical 
nodal  force  values.  Large  displacement  theory  was  used  in  the  nonlinear 
analyses. 
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RESULTS 


Experimental  Results 
Conventional  Gages 

Shear  moduli  as  determined  by  back-to-back  conventional  gages, 
back-to-back  full  section  gages  and  moir6/full  section  gages  are  listed  in 
Table  1.  Shear  modulus  was  calculated  by  linear  least  squares 
regression  of  data  between  1000  to  3000  (nominal)  microstrain. 

Statistics  are  reported  only  on  the  modulus  calculated  using  conventional 
gage  results  due  to  the  limited  number  of  specimens  evaluated  with  the 
other  strain  measurement  methods.  A  Student's  t  analysis  [27]  was  used 
to  compensate  for  the  small  sample  size. 

Ultimate  shear  strength  is  also  reported  in  Table  1.  Mean  strengths 
are  based  on  an  average  of  specimens  tested  at  NSWC  using  conventional 
gages  and  back-to-back  full  section  gages.  The  single  moird/full  section 
specimens  tested  at  VPI&SU  were  not  included  in  this  average. 

Recordings  of  the  U  and  V  displacement  fields  in  the  test  section  of 
each  of  the  four  moir6  specimens  at  specific  applied  average  shear  stress 
levels  are  shown  in  Figures  8  -11.  The  purity  of  shear  in  the  test  section 
is  characterized  by  the  lack  of  horizontal  fringe  gradients  in  the  V  field 
(Cy)  and  lack  of  vertical  fringe  gradients  in  the  U  field  (cx)-  Shear  strain 
fields  were  relatively  pure  away  from  the  notches  for  the  graphite  and 
glass  crossplys  and  the  glass  unidirectional  G23  specimen.  Axial  (cx) 
strains  at  the  notch  roots  indicated  the  presence  of  in-plane  bending 
possibly  due  to  the  low  flexural  stiffness  of  these  specimens,  particularly 
the  unidirectional  specimens.  This  is  similar  to  findings  by  Morton  et  al. 
[10]  for  90^  (fiber  orientation  relative  to  the  longitudinal  axis  of  the 
specimen)  graphite/epoxy  in -plane  shear  V-notched  beam  specimens. 
The  graphite/epoxy  unidirectional  G23  specimen  exhibited  the  least  pure 
shear  behavior  with  significant  axial  strains  in  the  test  section.  The 
glass/epoxy  crossply  specimen  in  Figure  11  exhibited  some  initial  Ex  at 
the  notch  roots  and  flanks  outside  the  test  section.  There  was  virtually 
zero  normal  strain  in  the  test  section  as  the  load  increased,  until  after 
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cracking  of  the  notch  roots/flanks  occurred.  Then  significant  Cx  and 
localized  shear  between  plys  became  evident  across  the  entire  test 
section. 

The  shear  strain  distribution  along  the  test  section  of  each  specimen 
was  calculated  from  the  moird  data,  and  is  shown  in  Figures  12-15.  The 
shear  strain  at  each  of  the  data  points  along  the  test  section  of  the 
specimens  has  been  normalized  with  the  mean  value  of  the  shear  strain 
for  the  moird  (front)  side  of  those  specimens.  Distributions  for  ail 
specimens  were  nonuniform,  however,  the  crossplys  exhibited  more 
uniform  behavior,  away  from  the  notch  roots,  than  the  unidirectional 
specimens. 

Moir6/Full  Section  and  Back-to-Back  Full  Section  Gages 

The  average  shear  strain  from  the  front  (moir6)  and  the  back  (full 
section  gage)  of  each  specimen  was  averaged  to  determine  the  true 
average  shear  strain  at  several  load  levels.  Plots  for  comparison  of  the 
moire/full  section  data,  back-to-back  full  section  gages,  and  back-to- 
back  conventional  gages  for  unidirectional  and  crossply  specimens  of 
both  material  systems  are  shown  in  Figures  16  -17.  Strain  differences 
between  the  front  and  back  of  the  specimen  were  observed  at  low  load 
levels,  but  decreased  with  increasing  load.  This  is  consistent  with 
observations  by  Morton  et  al.  (101  for  0°  in-plane  V-notched  beam 
specimens  which  showed  little  sensitivity  to  twisting  effects.  For  the 
unidirectional  G23  specimens,  the  moire/full  section  data  and  the  back- 
to-back  full  section  gage  data  are  in  excellent  agreement.  A  localized 
lower  stiffness  on  the  order  of  10-13%  is  calculated  with  the 
conventional  gage  data  which  does  not  account  for  the  nonuniform  strain 
distribution  across  the  entire  test  section.  A  similar  response  is  also 
apparent  in  the  crossply  graphite/epoxy  specimens,  especially  beyond 
0.7%  strain.  This  effect  is  not  as  severe  for  the  glass/epoxy  crossply 
specimens,  as  indicated  by  the  distributions  (Figure  15)  and  the 
relatively  good  agreement  between  the  moir6/full  section  data  and  the 
conventional  gage  data. 
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SEM 

Photomicrographs  of  graphite  and  glass  crossply  failures  are  shown 
in  Figure  18  (a)  and  (b)  respectively.  Shear  cracks  are  apparent 
between  0®  and  90°  laminae  and  tensile  cracks  are  apparent  within  the 
laminae.  The  AS4  specimens  also  showed  some  shear  cracks  within 
the  90®  laminae.  Photographs  of  characteristic  failures  of  the 
unidirectional  specimens  are  shown  in  Figure  19.  These  failures  arc 
reminiscent  of  those  observed  in  isotropic  specimens  as  reported  by 
Sullivan,  Kao  and  Van  Oene  [28]  and  consistent  with  assumptions  of 
transverse  isotropy  in  the  2-3  plane. 

Experimental/Analytical  Results 

Analysis  Versus  Conventional  Gaee  Results 

Analyses  with  linear  and  nonlinear  elastic  material  properties  are 
compared  to  data  from  the  062  back  to  back  strain  gages.  Gauss  point 
shear  strains  were  obtained  from  the  elements  physically  located 
beneath  the  strain  gage.  The  corresponding  shear  stress  was  the 
average  of  the  Gauss  point  shear  stresses  in  the  elements  from  notch 
root  to  notch  root  at  the  center  of  the  specimen.  Average  shear  stresses 
from  the  analysis  were  used  in  order  to  make  a  direct  comparison  with 
the  experimental  data  where  stresses  were  obtained  by  dividing  the 
load  by  the  cross-sectional  area  between  the  notches.  Boundary 
conditions  remained  the  same  for  analyses  with  both  linear  and 
nonlinear  elastic  material  properites.  No  attempt  was  made  to  model  the 
possibility  of  changing  boundary  conditions,  particularly  for  the  crossply 
specimens  with  measured  strains  near  3.5%. 

Figure  20  shows  predicted  through-thickness  shear  stress-strain 
curves  in  the  strain  gage  region  for  unidirectional  AS4/3501-6  and  S2 
glass/3501-6.  Only  slight  differences  are  observed  between  the  linear 
and  nonlinear  elastic  analyses.  Agreement  between  analysis  and  data  is 
excellent.  Corresponding  stress-strain  curves  for  the  cross-ply  laminates 
are  shown  in  Figure  21.  Significant  differences  are  observed  between 
the  linear  and  nonlinear  responses  beyond  1.5%  strain  for  AS4/3501-6 
and  1.0%  strain  for  S2  glass/3501-6.  The  measured  strains  exceed  those 
obtained  from  the  analysis  for  S2  giass/3501-6  beyond  1.0%  strain. 


Analysis  Versus  Moire/Full  Section  Gage  Data 

Comparisons  to  the  moird/full  section  gage  results  were  made  using 
both  the  average  stresses  and  strains  at  the  Gauss  points  of  the  elements 
from  notch  root  to  notch  root  in  the  specimen  center.  Analysis  results 
using  nonlinear  elastic  material  properties  agreed  well  with  the 
available  moird/full  section  and  back-to-back  full  section  data  points  for 
unidirectional  and  crossply  arrangements  of  both  materials  as  shown  in 
Figures  22  and  23.  Additional  moird  points  in  the  nonlinear  range  were 
not  obtained  for  the  crossply  laminates  due  to  difficulty  in  fringe 
resolution  at  higher  stress  levels. 

CONCLUSIONS 

Through-thickness  shear  V-notched  beam  shear  specimens  have  been 
investigated  both  experimentally  and  numerically  to  assess  strain  field 
uniformity  and  purity,  and  the  effect  of  3-D  nonlinear  vs.  3-D  linear  elastic 
material  properites  on  finite  element  analysis  results. 

Crossply  specimens  of  both  AS4  and  S2  glass  exhibited  uniform  strain 
distributions  away  from  the  notches  resulting  in  good  correlation  of 
conventional  gages  and  moir^/full  section  gage  data  at  low  stress  levels. 
Strain  distributions  in  the  AS4  and  S2  glass  unidirectional  specimens 
peaked  in  the  center  of  the  specimen  between  the  notches  so  that  strains 
measured  by  conventional  gages  exceeded  the  average  strain  across  the 
test  section.  This  resulted  in  a  lower  modulus  prediction  with  conventional 
gage  data  as  compared  to  moir^/full  section  data.  Significant  axial  strains 
existed  in  the  test  section  of  the  unidirectional  AS4  specimens,  while  the 
unidirectional  S2  glass  and  crossply  specimens  of  both  materials  showed 
fairly  pure  states  of  shear  in  the  test  section.  No  significant  strain 
differences  were  observed  between  the  front  and  back  of  the  specimen. 

Finite  element  analysis  results  using  a  3-D  nonlinear  elastic  material 
model  correlated  well  with  observed  responses  for  all  strain  measurement 
methods.  The  use  of  linear  elastic  material  properties  resulted  in  a  stiffer 
predicted  response  for  the  crossply  specimens  beyond  1%  and  1.5%  strain 
for  S2  glass  and  AS4  respectively.  The  inclusion  of  nonlinear  constitutive 
effects  are  recommended  in  the  analysis  of  these  specimens. 
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Ortsntatlon 


Material 


Strain  Maaauramant  Number  of 
Method  Spedmena 


Modulus  GPa(Msi)^  I  Strength  ^a(Msl)2 


IMdbadtanal 

AS4/3S01-6 

Conwanttonal 

3 

Z81  (0.40) 

1.9K 

GsgSpacknana 

MobalFid 

1 

3.14(0.45) 

32.45(4.69 

72% 

Baek-t»aack  Fid 

2 

3.14  (0.45) 

S2  glaBal3S01-6 

Corarartional 

3 

(0.00) 

13.0% 

30.75  (5.60) 

Moka/Fid 

1 

6.00(1.00) 

13.7% 

Baek-lo-Baefc  Fid 

2 

7.05(1.01) 

CroMpiy 

Sp#civvMns 

AS4/3S01>6 

Comraraianal 

4 

3.01  (0.56) 

2.3% 

01.37(13.08) 

MoiralFid 

t 

3.01  (0.56) 

1.0% 

8ack-lo43ack  Fid 

2 

3.01  (0.86) 

S2glaBar3S01-6 

Convanltonal 

4 

622(0.60) 

3.8% 

00.80(13.01) 

MokalFid 

t 

524  (0.75)3 

4.0% 

Baok-toSack  Fid 

2 

6.36(0.01) 

NOTES; 

1.  SMiwlci  oocnpultd  only  lor  aptcfciww  with  conwrticna)  oigw.  Nundm  In  par  ow<  (f or  both  modite  wd  Mrwi^N 
McNoeeeuraey  wtharaWaByof  9e%l»MdonaSiudwt  t  analyae  to  oornpanaala  ter  tie  amal  aamda  aiza  (27). 

2.  Moan  atranoih  k  taMd  on  an  awwao*  of  ttiangrte  oblainsd  from  ipadmana  wlh  oemartional  gagaa  and  baek-to-beek 
ful  aacHon  gaoat.  Al  ihaaa  apacfmani  wara  laatad  in  iha  aama  ffactura  at  NSWC. 

3.  Only  taro  mofrapoMa  at  low  toad  Itvatiwariavaadili  for  mod^caleidalionB.  Typtoaly.  aa  It*  toad  aea  fetcraaaad, 
thadacrapartcybalaraan  Irani  (mako)  and  back  (ltd  aaciton  gagaa)  daoraaaad.  Modubabaaad  on  ltd  aaelton  gagaa 
atona  to  6.45  QPa  (0.S2  MaQ. 


Table  1  Out-of-plane  shear  modulus  and  strength  values  for  unidirectional  and 
crossply  AS4/3501-6  and  S2  glass/3501-6  laminates 
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PROPERTY 

FJI 

BER  _ 

AS4 

S2GLASS 

Ei  (GPa) 

224 

88 

E2  (GPa) 

14 

88 

E3  (GPa) 

14 

88 

V12 

JZ 

.22 

''13 

.2 

.22 

''23 

.25 

.22 

Gi2(GPa) 

14 

35 

Gi3(GPa) 

14 

35 

G23(GPa) 

5.6 

35 

Subscripts  indicate  direction 
1  Along  fiber 
2.3  Transverse  to  fiber 


Table  2  AS4  graphite  and  S2  glass  fiber  properties  used  in  analysis 
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Property 

Epoxy  Resin 

3501-6 

Room  Temperature,  0%  Moisture 

K  (GPa) 

2.0 

G  (GPa) 

1.9 

P 

1.0 

N 

xq  (MPa) 

106 

Bulk  Modulus 
K-  E/(9-3E/G) 

Ramberg  Osgood  Shear  Model 

7  «  t/G  +  p{T(/G)(t/xo)^ 

t  Shear  stress 

7  Shear  strain 

G  initial  tangent  shear  modulus 

P.N.xq  Curvature  parameters 


Ted3le  3  3501-6  epoxy  matrix  properties  used  in  simplified  unit  ceil  model 
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ParaDOlfC  Mirror 


Figure  1  Scherratic  of  moire  setup  after  Guo  and  Ifju  [15] 
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(a) 


(b) 


Figure  3  Strain  gages  used  in  study  (a)  conventional  0/90 

stacked  gage,  (b)  full  section  gages  in  stacked  and 
side-by-side  configurations 
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Applied 

Load 


Wedges 


- Base 

■  ■'  —  ■  - - Plate 

12.70  cm  r 


Rgure  4  V-notched  beam  test  specimen  mounted  in  test  fixture  modeled 
after  fixture  developed  by  Adams  and  Walrath  [9] 
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3 


1 

7  «  Fiber  Direction 

I- -  1 - H 

Unit  Cell  (V  *  Fiber  Volume  Fraction) 


m  «  matrix 

For  component  1 1  components 

22.  12.  33.  23.  13 


Rgure  5  Simplified  unit  cell  model  based  on  spring  analogy  showing 
the  load  path  variations  between  fiber  and  matrix  dominant 
components  [24] 


T 

Figure  6  Predicted  through-thickness  shear  behavior  for 

unidirectional  (Q^)  AS4/3501-6  and  S2  glass/3501-6 
using  WY02D  and  the  simplified  unit  cell  model 


T 

Figure  7  Predicted  through-thickness  shear  behavior  for 

AS4/3501-6  and  S2  glass/3501-6  crossply  laminates 
using  the  simplified  unit  cell  model 
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Figure  8  U  (a)  and  V  (b)  moire  field  fringe  patterns  for  unidirectional  (G23) 

AS4/3501-6  specimens  loaded  at  23  MPa  and  7.7  MPa  respectively 


(a) 


(b) 


Figure  9  U  (a)  and  V  (b)  moir6  field  fringe  patterns  for  AS4/3501  -6  o’ossply 
specimen  at  a  load  level  of  14.3  MPa 
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Rgure  10  U  (a)  and  V  (b)  moirfi  field  fringe  patterns  for  unidirectional  (623) 
S2  glass/3501 -€  specimens  load^  at  21  MPa 


Figure  11  U  (a)  and  V  (b)  moir6  field  fringe  patttems  for  S2  glass/3501 -€ 

crossply  specimens  loaded  at  21 .8  MPa  and  1 1 .3  MPa  respectively 
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♦  t^3(23MPa) 


-  1  -0.5  0  0.5  1 


^max 

Figure  12  Normalized  shear  strain  distribution  along  the  test  section  of  a 
unidirectional  (G23)  AS4/3501-6  V-notched  beam  specimen 
measured  by  high  sensitivity  moir6  at  two  load  levels 


Rgure  13  Normalized  shear  strain  distribution  along  the  test  section  of  a 
unidirectional  (G23)  S2  glass/3501-6  V-notched  beam  specimen 
measured  by  high  sensitivity  moir4  at  several  load  levels 


30 


-1  -0.5  0  0.5  1 

^max 

Figure  14  Normalized  shear  strain  distribution  along  the  test  section  of  a 
crossply  AS4/3501-6  V-notched  beam  specimen  measured  by 
high  sensitivity  moir6  at  two  load  levels 


*^max 

Rgure  15  Normalized  shear  strain  distribution  along  the  test  section  of  a 
crossply  S2  glass/3501  <6  V-notched  beam  specimen  measured 
by  high  sensitivity  rnoir^  at  two  load  levels 
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Shear  Stress,  MPa 


Figure  16  Shear  stress-strain  response  for  AS4/3501-6  and 

S2  glass/3501-6  unidirectional  (G23)  specimens  instrumented 
with  conventional  gages  vs.  moir6  and  full  section  gage  data 


Rgure  17  Shear  stress-strain  rwponse  for  AS4/3501-6  and 

S2  glass/3501-6  crossply  specimens  instrumented  with 
conventional  gages  v$.  moir6  and  full  section  gage  data 
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(a)  AS4/3501-6 


Figure  18  Photomicrographs  of  AS4/3501-6  and  S2  glass/3501-6 
crossply  specimen  failures 


Shear  Stress.  MPa  Shear  Stress.  MPa 


Figure  20  Shear  stress-strain  response  for  unidirectionai  (G23)  AS4/3501-6 
and  S2  glass/3501-6  specimens  instrumented  with  conventional 
gages  (062)  vs.  predicted  values  based  on  finite  element 
analyses  with  linear  and  nonlinear  elastic  material  prcoerties 


Rgure  21  Shear  stress-strain  response  of  crossply  specimens  predicted 
by  finite  element  analyses  with  linear  and  nonlinear  elastic 
material  prooerties.  and  measured  by  conventional  (062)  gages 
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Figure  22  Average  shear  stress-strain  response  for  AS4/3501  -6  and 

S2  gfassA350l-6  unidirectional  (G23)  spedmens  measured  with 
moir6  interferometry  and  full  section  gages  vs.  predicted  values 
based  on  finite  element  analyses  with  nonlinear  elastic  material 
properties 


Rgure  23  Average  shear  stress-strain  response  for  AS4/3501-6  and 
S2  glass/3501'6  crossply  specimens  measured  with  moir6 
interferometry  and  full  section  gages  vs.  predicted  vdues 
based  on  finite  element  analyses  with  nonlinear  elastic  material 
properties 
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